Molecular markers for variation in spawning date in a hatchery population of rainbow trout (Oncorhynchus mykiss).

We examined the distribution of alleles at 63 microsatellite loci distributed across 29 linkage groups in broodstock females from a commercial population of rainbow trout spawning on different dates throughout the season (August to January). A total of 368 females, 184 and 117 females from each of the tail-ends of the spawning distribution and a subsample of 67 females spawning in the middle, were used to detect marker-trait associations. Twenty-one loci in a subset of genomic regions (RT-5, 7, 8, 10, 12, 14, 15, 22, 23, 24, 25, 29, 30, and 31) were significantly associated with variation in spawning date. Many of these markers localize to regions with known spawning date quantitative trait loci based on previous studies. An individual assignment analysis was used to test how well the molecular data could be used to assign individuals to their correct spawning group, and markers were given a ranking reflecting their contribution to the accuracy of assignment. The top 15 ranked markers were successful at assigning the majority of females to the correct spawning group based on genotype with an average accuracy of 76 %. The most likely genes that could contribute to these differences in spawning date are discussed. Together, these data indicate that the loci could be incorporated into a selection index with phenotype data to increase the accuracy of selection for spawning date.

Comparative analysis of genetic parameters and quantitative trait loci for growth traits in Fraser strain Arctic charr (Salvelinus alpinus) reared in freshwater and brackish water environments.

To determine the potential for genetic improvement in Fraser strain Arctic charr (AC, Salvelinus alpinus), we calculated genetic parameters for BW and condition factor (K) and tested if previously identified QTL for these traits were detectable across a commercial broodstock reared in both freshwater (FRW) and brackish water (BRW). Individuals from 30 full-sib families were reared up to 29 mo of age in FRW and BRW tanks at a commercial facility. Heritability for BW and K was moderate in FRW (0.29 to 0.38) but lower in BRW (0.14 to 0.17). Genetic correlations for BW across environments were positive and moderate (0.33 to 0.67); however, equivalent K correlations were very weak (0.24 to 0.37). We identified a single BW QTL with experimentwide effects on linkage group AC-8, 4 BW QTL (AC-4, -13, -14, and -19), and 3 K QTL (AC-4, -5, and -20) with chromosomewide effects across families. Notably, the QTL on AC-8 had significant effects with BW at 3 out of 4 sampling dates in FRW and had significant allelic phase disequilibrium with BW across families, suggesting a tight coupling of the marker region to the QTL in this population. Body weight QTL were identified on AC-4 in both FRW and BRW environments and AC-4 was the only linkage group with a detectable QTL for both K and BW. Modest consistency of some QTL effects as well as moderate heritability in both environments suggests that there is some potential for genetic improvement of growth in this species even though gene × environment interactions are high.

Whole genome duplication: challenges and considerations associated with sequence orthology assignment in Salmoninae.

To illustrate some of the challenges and considerations in assigning correct orthology necessary for any comparative genomic investigation among salmonids, sequence data from the non-coding regions of different chromosomes in three members of the subfamily Salmoninae, rainbow trout Oncorhynchus mykiss, Atlantic salmon Salmo salar and Arctic charr Salvelinus alpinus, were compared. By analysing c. 55 distinct loci, corresponding to c. 142 kbp sequence information per species, 18 duplicated patterns representative of the two sequential rounds of teleost-specific whole genome duplications (i.e. 3R and 4R WGD) were identified. Sequence similarities between the 4R paralogues were c. 90%, which was slightly lower than those of the 4R orthologues and c. 60% for the 3R products. Through careful examination of the sequence data, however, only 14 loci could reliably be assigned as true orthologues. Locus-specific trees were constructed through maximum parsimony, maximum likelihood and neighbour-joining methods and were rooted using the information from a close relative, lake whitefish Coregonus clupeaformis. All approaches generated congruent trees supporting the {Coregonus [Salmo (Oncorhynchus, Salvelinus)]} topology. The general phenotypic characteristics of sequences, however, were highly suggestive of the basal position of Oncorhynchus, raising the hypothesis of an accelerated rate of nucleotide evolution in this species.

Phenotypic plasticity, heterochrony and ontogenetic repatterning during juvenile development of divergent Arctic charr (Salvelinus alpinus).

Phenotypic plasticity is a developmental process that plays a role as a source of variation for evolution. Models of adaptive divergence make the prediction that increasing ecological specialization should be associated with lower levels of plasticity. We tested for differences in the magnitude, rate and trajectory of morphological plasticity in two lake populations of Arctic charr (Salvelinus alpinus) that exhibited variation in the degree of resource polymorphism. We reared offspring on diet treatments that mimicked benthic and pelagic prey. Offspring from the more divergent population had lower levels of morphological plasticity. Allometry influenced the rate of shape change over ontogeny, with differences in rate among ecomorphs being minimal when allometric variation was removed. However, plasticity in the spatial trajectory of development was extensive across ecomorphs, both with and without the inclusion of allometric variation, suggesting that different aspects of shape development can evolve independently.

The genetic architecture of embryonic developmental rate and genetic covariation with age at maturation in rainbow trout Oncorhynchus mykiss.

The genetic architecture underlying variation in embryonic developmental rate (DR) and genetic covariation with age of maturation (MAT) was investigated in rainbow trout Oncorhynchus mykiss. Highly significant additive parental effects and more limited evidence of epistatic effects on progeny hatching time were detected in three diallel sets of families. Genome scans with an average of 142 microsatellite loci from all 29 linkage groups in two families detected significant quantitative trait loci (QTL) for developmental rate on RT-8 and RT-30 with genome-wide and chromosome-wide effects, respectively. The QTL on linkage group RT-8 explained 23·7% of the phenotypic variation and supports results from previous studies. The co-localization of QTL for both DR and MAT to several linkage groups and the observation that alleles associated with faster developmental rate were found significantly more often in early maturing rather than typical and later maturing male ancestors supports the hypothesis of genetic covariation between DR and MAT. The maturation background and schedule of additional sires, however, did not have a consistent association with their progeny hatching times, suggesting that other genetic, environmental and physiological effects contribute to variation in these life-history traits.

Evolution of adaptive diversity and genetic connectivity in Arctic charr (Salvelinus alpinus) in Iceland.

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.

A framework for power and sensitivity analyses for quantitative genetic studies of natural populations, and case studies in Soay sheep (Ovis aries).

Studies of the quantitative genetics of natural populations have contributed greatly to evolutionary biology in recent years. However, while pedigree data required are often uncertain (i.e. incomplete and partly erroneous) and limited, means to evaluate the effects of such uncertainties have not been developed. We have therefore developed a general framework for power and sensitivity analyses of such studies. We propose that researchers first generate a set of pedigree data that they wish to use in a quantitative genetic study, as well as data regarding errors that occur in that pedigree. This pedigree is then permuted using the data regarding errors to generate hypothetical 'true' and 'assumed' pedigrees that differ so as to mimic pedigree errors that might occur in the study system under consideration. Phenotypic data are then simulated across the true pedigree (according to user-defined genetic and environmental covariance structures), before being analysed with standard quantitative genetic techniques in conjunction with the 'assumed' pedigree data. To illustrate this approach, we conducted power and sensitivity analyses in a well-known study of Soay sheep (Ovis aries). We found that, although the estimation of simple genetic (co)variance structures is fairly robust to pedigree errors, some potentially serious biases were detected under more complex scenarios involving maternal effects. Power analyses also showed that this study system provides high power to detect heritabilities as low as about 0.09. Given this range of results, we suggest that such power and sensitivity analyses could greatly complement empirical studies, and we provide the computer program PEDANTICS to aid in their application.

Genomic organization of the IGF1, IGF2, MYF5, MYF6 and GRF/PACAP genes across Salmoninae genera.

Whole-genome duplication in the ancient ray-finned fish and subsequent tetraploidization in the ancestor to the salmonids have complicated genomic and candidate gene studies in these organisms as many genes with multiple copies are present throughout their genomes. In an attempt to identify genes with a potential influence on growth and development, we investigated the genomic positions of insulin-like growth factors 1 and 2 (IGF1, IGF2), myogenic factors 5 and 6 (MYF5, MYF6) and growth hormone-releasing factor/pituitary adenylate cyclase-activating polypeptide (GRF/PACAP) in three salmonid species: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar) and Arctic charr (Salvelinus alpinus). Our results suggest a tight association between the IGF1, MYF5 and MYF6 genes in all three species. We further localized the duplicated copies of IGF1 to the homeologous linkage groups RT-7/15 in rainbow trout and AC-3/24 in Arctic charr, and the two copies of MYF6 to homeologous linkage groups AS-22/24 in Atlantic salmon. Localization of GRF/PACAP to RT-7, AS-31 and AC-27 and IGF2 to RT-27, AS-2 and AC-4 in rainbow trout, Atlantic salmon and Arctic charr respectively is consistent with previously reported homologies among these chromosomal segments identified using other genetic markers. However, localization of the second copy of GRF/PACAP to RT-19 and AC-14 and the duplicated copy of IGF2 to AC-19 suggest a possible new homology/homeology between these chromosomes. These results might also be an indication of a more ancient polyploidization event that occurred deep in the ray-finned fish lineage.

Genome organization of glutamine synthetase genes in rainbow trout (Oncorhynchus mykiss).

Unlike mammals, bony fish appear to possess multiple genes encoding glutamine synthetase (GS), the nitrogen metabolism enzyme responsible for the conversion of glutamate and ammonia into glutamine at the expense of ATP. This study reports on the development of genetic markers for each of the four isoforms identified thus far in rainbow trout (Oncorhynchus mykiss) and their genome localization by linkage mapping. We found that genes coding for GS01, GS02, GS03, and GS04 map to four different linkage groups in the trout genome, namely RT-24, RT-23, RT-08, and RT-13, respectively. Linkage groups RT-23 and RT-13 appear to represent distinct chromosomes sharing duplicated marker regions, which lends further support to the previous suggestion that GS02 and GS04 may be duplicate gene copies that evolved from a whole-genome duplication in the trout ancestor. In contrast, there is at present no further evidence that RT-24 and RT-08 share ancestrally homologous segments and additional genomic studies will be needed to clarify the evolutionary origin of genes coding for GS01 and GS03.

Multiple hyperechoic testicular lesions are a common finding on ultrasound in Cowden disease and represent lipomatosis of the testis.

Cowden disease (CD) is a genetic disease associated with multiple hamartomas and malignant neoplasms. During investigations for possible subnormal fertility, a series of eight males with CD underwent ultrasound scanning of their testes. Our findings detail the seven adult patients that were found to have multiple and bilateral testicular lesions on ultrasound. These lesions in Cowden's patients represent a newly described testicular pathology - lipomatosis of the testis. Here we detail the radiological findings. Ultrasound findings showed multiple (estimated 40+) discreet lesions randomly scattered thoughout the parenchyma of each affected testis. The lesions themselves had heterogeneous echo texture, but all being hyperechoic and with a variable size from 1 mm to 6 mm. One patient underwent MR examination of the lesions in which a high T(1) signal was seen, but no useful characteristic features were identifiable.

The candidate gene, Clock, localizes to a strong spawning time quantitative trait locus region in rainbow trout.

We applied a candidate gene mapping approach to an existing quantitative trait loci (QTL) data set for spawning date in rainbow trout (Oncorynchus mykiss) to ascertain whether these genes could potentially account for any observed QTL effects. Several genes were chosen for their known or suspected roles in reproduction, circadian, or circannual timing, including salmon-type gonadotropin-releasing hormone 3A and 3B (GnRH3A and GnRH3B), Clock, Period1, and arylalkylamine N-acetlytransferase-1 and -2 (AANAT-1 and AANAT-2). Genes were sequenced, and polymorphisms were identified in parents of two rainbow trout mapping families, one of which was used previously to detect spawn timing QTL. Interval mapping was used to identify associations between genetic markers and spawning date effects. Using a genetic map that was updated with 574 genetic markers (775 total), we found evidence for 11 significant or suggestive QTL regions. Most QTL were only localized within one of the parents; however, a strong QTL region was identified in both female and male parents on linkage group RT-8 that explained 20% and 50% of trait variance, respectively. The Clock gene mapped to this region. Period1 mapped to a region in the female parent associated with a marginal effect (P = .056) on spawn timing. Other candidate genes were not associated with significant QTL effects.

Aplasia of the parotid gland in Down syndrome.

Salivary gland aplasia has not to our knowledge been previously reported in association with Down syndrome. We present a case of bilateral parotid aplasia in a patient with Down syndrome. Clinically he had aplasia of the major salivary glands and symptoms of xerostomia. Thirteen other family members over three generations were examined, and all had functional parotid glands. We reviewed publications about Down syndrome and salivary aplasia, together with the data regarding his other clinical problems and family background. His oral problems were inadequate plaque control, dental caries, and erosion of the teeth.

Sex-linked quantitative trait loci for thermotolerance and length in the rainbow trout.

We hypothesized that correlation between growth traits and upper thermal tolerance (UTT) in rainbow trout (Oncorhynchus mykiss) might be explained by quantitative trait loci (QTL) localized to the same linkage groups. Microsatellites on three autosomal linkage groups carrying UTT QTL in rainbow trout were tested for associations with fork length (FL) and condition factor (K) in half-sib families of outbred rainbow trout and in backcrosses of trout lines selected on UTT. Additionally, we used a sex-linked microsatellite (OmyFGT19TUF) to test for marker-trait associations at the sex chromosomes. The sex-linked marker OmyFGT19TUF was significantly associated with FL and UTT, accounting for up to 9.6% and 9.7% of variance in these traits, respectively. Male advantages in FL (and, to a lesser extent, UTT) relative to their female sibs were dependent on the origin of the Y chromosome and thus varied among grandsire lines. However, males had higher K in a manner unrelated to Y chromosomal origin, suggesting a partially sex-limited expression of this trait. Omy325UoG was significantly associated with K in one of the outbred half-sib families, but no other significant autosomal marker-trait associations were detected. Our findings illustrate minor evidence that correlation between UTT and FL is partially determined by one or more sex-chromosomal QTL.

QTL for body weight and condition factor in Atlantic salmon (Salmo salar): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Arctic charr (Salvelinus alpinus).

Genotypes at 91 microsatellite loci in three full-sib families were used to search for QTL affecting body weight (BW) and condition factor in North American Atlantic salmon (Salmo salar). More than one informative marker was identified on 16-18 linkage groups in each family, allowing at least one chromosomal interval to be analyzed per linkage group. Two significant QTL for BW on linkage groups AS-8 and AS-11, and four significant QTL for condition factor on linkage groups AS-2, AS-5, AS-11, and AS-14 were identified. QTL for both BW and condition factor were located on linkage groups AS-1, 6, 8, 11, and 14 when considering both significant and suggestive QTL effects. The largest QTL effects for BW (AS-8) and for condition factor (AS-14) accounted for 20.1 and 24.9% of the trait variation, respectively. Three of the QTL for BW occur on linkage groups where similar effects have been detected on the homologous regions in either rainbow trout (Oncorhynchus mykiss) or Arctic charr (Salvelinus alpinus).

Linkage arrangement of Na,K-ATPase genes in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss).

As part of our efforts to characterize Na,K-ATPase isoforms in salmonid fish, we investigated the linkage arrangement of genes coding for the alpha and beta-subunits of the enzyme complex in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss). Genetic markers were developed from four of five previously characterized alpha-subunit isoforms (alpha1b, alpha1c, alpha2 and alpha3) and four expressed sequence tags derived from yet undescribed beta-subunit isoforms (beta1a, beta1b, beta3a and beta3b). Sex-specific linkage analysis of polymorphic loci in a reference meiotic panel revealed that Na,K-ATPase genes are generally dispersed throughout the rainbow trout genome. A notable exception was the colocalization of two alpha-subunit genes and one beta-subunit gene on linkage group RT-12, which may thus share a conserved orthologous segment with linkage group 1 in zebrafish (Danio rerio). Consistent with previously reported homeologous relationships among the chromosomes of the rainbow trout, primers designed from the alpha3-isoform detected a pair of duplicated genes on linkage groups RT-27 and RT-31. Similarly, the evolutionary conservation of homeologous regions on linkage groups RT-12 and RT-16 was further supported by the map localization of gene duplicates for the beta1b isoform. The detection of homeologs within each gene family also raises the possibility that novel isoforms may be discovered as functional duplicates.

Population genetic structure of Arctic charr, Salvelinus alpinus from northwest Europe on large and small spatial scales.

To examine the population genetic structure of lake-resident Arctic charr, Salvelinus alpinus from northwest Europe on multiple spatial scales, 2367 individuals from 43 lakes located in three geographical regions (Iceland, the British Isles and Scandinavia) were genotyped at six microsatellite loci. On a large scale, data provided little evidence to support clustering of populations according to geographical region. Hierarchical analysis of molecular variance indicated that, although statistically significant, only 2.17% of the variance in allelic frequencies was partitioned at the among-region level. Within regions, high levels of genetic differentiation were typically found between lakes regardless of the geographical distance separating them. These results are consistent with the hypothesis of rapid postglacial recolonization of all of northwest Europe from a single charr lineage, with subsequent restriction of gene flow. On a smaller scale, there was evidence for close genetic relationships among lakes from within common drainage basins in Scotland. Thus, interlake genetic structure reflects localized patterns of recent (or contemporary) gene flow superimposed onto a larger scale structure that is largely a result of historical processes. There was also evidence for widespread genetic structuring at the within-lake level, with sympatric populations detected in 10 lakes, and multilocus heterozygote deficits found in 23 lakes. This evidence of the Wahlund effect was found in all lakes known to contain discrete phenotypic morphs, as well as many others, suggesting that morphs may often represent separate breeding populations, and also that the phenomenon of polymorphism in this species may be more widespread than is currently realized.

A genetic linkage map for Arctic char (Salvelinus alpinus): evidence for higher recombination rates and segregation distortion in hybrid versus pure strain mapping parents.

We constructed a genetic linkage map for Arctic char (Salvelinus alpinus) using two backcrosses between genetically divergent strains. Forty-six linkage groups (expected = 39-41) and 19 homeologous affinities (expected = 25) were identified using 184 microsatellites, 129 amplified fragment length polymorphisms (AFLPs), 13 type I gene markers, and one phenotypic marker, SEX. Twenty-six markers remain unlinked. Female map distance (9.92 Morgans) was substantially higher than male map distance (3.90 Morgans) based on the most complete parental information (i.e., the F1 hybrids). Female recombination rates were often significantly higher than those of males across all pairwise comparisons within homologous chromosomal segments (average female to male ratios within families was 1.69:1). The female hybrid parent had significantly higher recombination rates than the pure strain female parent. Segregation distortion was detected in four linkage groups (4, 8, 13, 20) for both families. In family 3, only the largest fish were sampled for genotyping, suggesting that segregation distortion may represent regions possessing influences on growth. In family 2, almost all cases showing segregation distortion involved markers in the female hybrid parent.

Comparison of GNRH3 genes across salmonid genera.

Genomic sequences of gonadotropin-releasing hormone genes were amplified and examined for sequence divergence among members of three different genera of the subfamily Salmoninae: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus). Sequences of GNRH3A and GNRH3B (formerly known as sGnRH1 and sGnRH2) were 97-99% similar in coding regions and 94-98% similar in non-coding regions among genera, but comparisons within species between GNRH3A and GNRH3B were only 90-92% similar in coding regions and 83-89% similar in non-coding regions. Polymorphisms in the parents of mapping families for each species allowed for linkage mapping of the GNRH3B gene in all three species and the GNRH3A gene in rainbow trout. GNRH3B maps to linkage group 6 in rainbow trout, linkage group 16 in Atlantic salmon and linkage group 25 in Arctic charr. GNRH3A mapped to linkage group 30 in rainbow trout.

Selective and genetic constraints on the evolution of body size in a stream-dwelling salmonid fish.

To examine constraints on evolution of larger body size in two stunted populations of brook charr (Salvelinus fontinalis) from a single river in Cape Race, Newfoundland, Canada, we measured viability selection acting on length-at-age traits, and estimated quantitative genetic parameters in situ (following reconstruction of pedigree information from microsatellite data). Furthermore we tested for phenotypic differentiation between the populations, and for association of high juvenile growth with early maturity that is predicted by life history theory. Within each population, selection differentials and estimates of heritabilities for length-at-age traits suggested that evolution of larger size is prevented by both selective and genetic constraints. Between the populations, phenotypic differentiation was found in length-at-age and age of maturation traits, whereas early maturation was associated with increased juvenile growth (relative to adult growth) both within and between populations. The results suggest an adaptive plastic response in age of maturation to juvenile growth rates that have a largely environmental basis of determination.